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.Abstract 

■We study the sensitivity of future medium baseline reactor antineutrino experiments on the neutrino mass hierarchy. By using the 
Istandard^^ analysis, we find that the sensitivity depends strongly on the baseline length L and the energy resolution (6E/Ef' = 

where a and b parameterize the statistical and systematic uncertainties, respectively. The optimal length is 
■found to be L ~ 40 - 55 km, the larger resolution the shorter optimal L. The running time needed to determine the mass hierarchy 
^depends strongly on the energy resolution. The 3cr determination is possible after 10 years of running with a = 3% and b - 0.5%, 
■which can be reduced to 4 years if a = 2% and b = 0.5%, for a 5 kton detector (with 12% weight fraction of free proton) placed at 
'L ~ 50 km away from a 20 GWth reactor. This type of experiment can measure the relevant mixing parameters accurately, achieving 
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6 sin^ 26i|2 ~ 4 X 10"^ 6(ml - m\) ~ 0.03 x 10"^eV^ and 5|m^ - m^^ \ ~ 0.007 x lO^^eV^, after 5 years 
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'1. Introduction 

. Now that a large has been measured at Daya Bay [1] 
"and RENO [2] experiments accurately, neutrino physics en- 
.ters a new era. One of the next challenges is determination 
'of the mass hierarchy. Many ideas have been proposed, such 
.as long baseline accelerator based neutrino oscillation [3-5], 
■atmospheric neutrino [6], supernova neutrino [7], neutrino-less 
.double-beta decay [8] and medium baseline reactor antineutrino 
■experiments [9-15]. 

, Among them, the medium baseline reactor antineutrino ex- 
periment has stimulated various re-evaluations of its physics 
.potential and sensitivity recently. Some works utilize 
■the Fourier transform technique [16-18], first discussed in 
,Refs. [11-13], to distinguish the mass hierarchy. The main 
advantage of this technique is that the mass hierarchy can be 
determined without precise knowledge of reactor antineutrino's 
spectrum, the absolute value of the large mass-squared differ- 
ence lAm^jl, and the energy scale of a detector Although in- 
teresting and attractive, this technique is somewhat subtle to in- 
corporate the uncertainties of the mixing parameters and to es- 
timate its sensitivity to the mass hierarchy. On the other hand, 
some works adopt the;\f^ analysis [15, 18, 19] and new mea- 
sure based on Bayesian approach[20]. These methods utilize 
all available information from experiments, and it is straightfor- 
ward to incorporate the uncertainties to evaluate the sensitivity. 
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providing robust and complementary results to the Fourier tech- 
nique. 

In this paper, we analyze the sensitivity of medium base- 
line reactor antineutrino experiments to the mass hierarchy for 
the baseline length of 10-100 km and the energy resolution 

(SEjEf = [a/y/E/Mevf + b^ in the range 2% < a < 6% 
and b < 1% with the standard analysis. The optimal base- 
line length and the expected statistical uncertainties of neutrino 
parameters, sin^ 20i2, sin^ 200, Amjj and Am^j, are also esti- 
mated. 

2. Reactor antineutrino flux 

In this section, we briefly discuss the evaluation of how many 
electron antineutrinos Vg would be detected at a far detector with 
a medium baseline length from a reactor 

In a nuclear reactor, antineutrinos are mainly produced via 
beta decay of the fission products of the four radio-active iso- 
topes, 235f/, 238f; 239p^^ 241 jj^g f^^gji number 

of antineutrinos produced per fission depends on their energy 
Ey [22] 

4>(Ey) = /235f;exp (0.870 -0.160£v-0.091£^) 
-1- /2_«p_^ exp (0.896 - 0.239£v - 0.098 1£^) 
H- /23sy exp (0.976 - 0.162£v - 0.0790£^) 
H- /24,p_^ exp (0.793 - 0.080£v - 0.1085 E^) , (1) 



'Precisely speaking, there ai'e contributions from other isotopes such as 
"P,, and ^^^Pu, but their contributions are of the order of 0.1% or less [21]. 



where fk denotes the relative fission contribution of the isotope 



A: in a reactor fuel, derived from the fission rate AT/^"* {s ') of 



isotope k as 



N; 



■fiss 



rfiss * 



(2) 



Although fic varies over time as the fuel is burned, it can be ap- 
proximated for this type of experiments with the average value 
of the relative fission contributions: /235U = 0.58, f239p^ = 
0.30, /asy = 0.07 and /24ip„ = 0.05 [12]. The event rate of 
antineutrinos with energy Ey (MeV) at a reactor of Pi GWth) 
thermal power is then expressed as 



dN 
dEy 



fk^k 



4>(Ey)x 6.24 X 10 



21 



(3) 



where ek is the released energy per fission of the isotope k: 
e2,5u = 201.7 MeV, e2,9p^ = 210.0 MeV, eiisj; = 205.0 MeV 
and £241 = 212.4 MeV |23|. The numerical factor comes from 
unit conversion, 1 GW/MeV = 6.24 x 10^^ 

This rate is then modulated by oscillation. The Vg survival 
probabiUty is expressed as 

2 



(=1 



jCXp 



2Ei 



= 1 - cos'' 013 sin^ 2012 sin^ (A21) 

- cos^ 012 sin^ 2013 sin^ (A31) 

- sin^ 012 sin^ 20i3 sin^ (A32) , 



(4) 



where Uej is the neutrino mixing matrix element relating the 
electron neutrino to the mass eigenstate v,. The variables m; and 
Ei are the mass and energy of the corresponding mass eigen- 
state, while 6ij represent the neutrino mixing angles. The oscil- 
lation phase Aij are defined as, 

AmlL 



{Anijj = ntj - ntj). 



(5) 



with a baseline length L. We have neglected the matter effect 
because it is small enough for the energy range and the baseline 
lengths we concern in this study [24]. In obtaining the second 
line of (4) we have also ignored the tiny energy difference be- 
tween the three mass eigenstates, Ey ~ Ei ~ E2 ~ £3. 

To make the eflTect of mass hierarchy more clear, we would 
like to rewrite eq. (4) as. 

Pee = 1- cos^ 013 sin^ 2012 sin^ (A21) 
-sin2 20i3sin2(|A3i|) 
- sin^ 012 sin^ 20i3 sin^ (A21) cos (2IA31 1) 

+ sin2 2013 sin (2A21) sin (2IA31 1) , (6) 

where only the last term depends on the mass hierarchy, which 
takes the plus and minus sign, respectively, for normal (NH) 
and inverted hierarchy (IH), 



.3 - m\ > (NH) 
ml-m\<0 (IH). 



(7a) 
(7b) 



It is clear from eq. (6) that the survival probability is most sensi- 
tive to the mass hierarchy when | sin(2A2i)| = I, or equivalently 

2A2i=(2«-l)^ (n= 1,2,3,---), (8a) 

and has no sensitivity at 

2A2i=n;r (n = 0, 1, 2, 3, • • ■), (8b) 

where sin(2A2i) - 0. For example, at L = 50 km, the condition 
(8a) for n - \ and 2 is satisfied at ~ 6 MeV and 2 MeV, re- 
spectively. The last term in eq. (6) contributes with the opposite 
sign at these first and second maxima. In between, it vanishes 
and changes its sign at Ziy = 3 GeV, corresponding to « = 1 in 
(8b). It is this sign change that plays an important role for the 
mass hierarchy determination, which will be further discussed 
in the next section. 

Similar as the current reactor experiments, such as Daya 
Bay [1], RENO [2] and Double Chooz [25], future medium 
baseline reactor antineutrino experiments can also use free pro- 
tons as targets to detect electron antineutrinos via the inverse 
neutron beta decay (IBD) process. 



Ve +P 



+ n, 



(9) 



where p and n are the proton and the neutron, respectively. The 
threshold neutrino energy of this process is Efhr ~ m„-mp+me, 
and the cross section is [26], 



-42 2 

cm^. 



(10) 



where Eg and pe are the energy and momentum of the positron, 
neglecting the kinetic energy of the proton and the neutron for a 
MeV scale antineutrino. The positron's energy is roughly Ee ~ 
Ey - {nin - nip). 

The produced positron then interacts with scintillator, con- 
verting its kinetic energy to photons. Eventually, the positron 
annihilates with an electron in the detector and emits two 0.5 
MeV photons. The total energy of those photons are then ac- 
cumulated as the visible energy, £'v/j, which is a sum of the 
positron's total and one electron's rest energies. 



-Ee + Me ~ {Ey - Q.%)MqW . 



(11) 



Finite energy resolution of the detector then distorts the 
true visible energy, Eyi^, to the finally observed one, E"'^^. 
This effect can be modeled by a detector response function 
G(£,/, - 5E,i,) with the energy resolution SEyi^. In this 
study, we take the normahzed gaussian function as the response 
function, i.e.. 



G{Eyis - E°f \ SEyis) 



1 



^l2n6Ey 



■ exp < 



2((5£w.)2 



The detector energy resolution [14] 



6E,^ 



yfE~jMSV 



(12) 



(13) 



is composed of two parts. The first term in the square-root rep- 
resents the statistical uncertainty, and the second one gives the 
systematic uncertainty [27]. The observed antineutrino distri- 
bution by a detector with Np free protons after an exposure time 
T can then be expressed as 



dN N„T 



dE' 



obs 



NpT r°° dN 



X criBD(Ey) G(E, - O.SMeV - E^f' 6E,i,). (14) 



3. The sensitivity to the mass hierarchy 

After obtaining the energy distribution of reactor antineutri- 

nos, we would like to estimate the sensitivity of determining the 
mass hierarchy using the standard analysis [10, 13, 15, 18, 
19]. 

To set the stage, we introduce the;^'^ function as 



.2 

Stat' 



(15) 



The first term summarizes the prior knowledge on mixing pa- 
rameters. In reactor antineutrino experiments, these are the 
naixing angles, sin^ 20i2 and sin^ 26is, and the two mass-square 
differences, Am^^ and | Awjj j |, whose contributions look like, 

2 _ f (sin^ 2012)^" - (sin^ 2012) 

r (sin^ 26*13) - (sin^ 26li3)'"^"' 
1 (5sin^20i3 

^1 liKi^i j ■ 



(16) 



The input values Y'"P"' and their uncertainties SY are listed in 
Table 1. 

The reactor antineutrino flux, IBD cross section, fiducial vol- 
ume and weight fraction of free proton can all be combined into 
a single overall factor Consequently, their contributions to the 

function can be represented by a single term as. 



X] 



( ff'' 

J sys 



sys 



r input 
■ Jsys 



Sfsy 



(17) 



where f^f"' = 1, and Sf^ys = 0.03. 

The third term in (15) represents the statistical fluctuation. 
When we introduce bining w.r.t. E°f^\ it looks like 



Y 


sin^ 26n 


sin^ 2^13 


Am^i eV^ 




fsys 


yinput 

5Y 


0.857 
0.024 


0.089 
0.005 


7.50 X 10-5 
0.20 X 10-5 


2.32 X 10-3 
0.1 X 10-3 


1 

0.03 



Table 1: The input values Y'"'"" and their uncertainties 6Y taken from Refs. [1, 
28]. The uncertainty of sin^ 19^ can be 5% or less after 3 years running of 
Daya Bay experiment [29]. 



with the summation running over all the bins. Here, N^^'^'^ is 
the event number for the ;,/, bin when the hierarchy is NH (IH), 
while N^" is the theoretical prediction of the event number ei- 
ther with right or wrong mass hierarchy, calculated as a function 
of the four model parameters and the normalization factor /jv,, 
which are all varied under the constraints of (16) and (17). In 
this study we prepare the data by using eq.(14) with the 

input values of the five parameters for each mass hierarchy. 

In the limit of infinitely many events, the bin size can be re- 
duced to zero, and the sum (18) can be replaced by an integral. 



y2 

A-stc 
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dE' 



obs 
vis 



dN Y I dN \ 



NH(IH) 



^NH(IH) 



(19) 



where = 1.8 MeV and E,„ax - 8 MeV are the lower and 
upper limits of the observed energy used to evaluate the^^ func- 
tion, respectively. Although a finite bin size is required for ac- 
tual experiments, we adopt this zero-bin-size limit as measure 
of the maximum sensitivity. We then define Ax^ as 



Ax X Xmin ' 



(20) 



where xlii„ is flie minimum of which is obviously zero in 
our approximation of neglecting statistical fluctuations in data, 
AT™. When wrong mass hierarchy is assumed in the fit, 
the minimum of Ax^, (Ax^)min, will deviate from zero, and 
the wrong mass hierarchy can be rejected with significance 

y/(Ax^)min^- 

4. Results 

In this section, we discuss the sensitivity to the mass hier- 
archy, the optimal length and the statistical uncertainties of the 
neutrino parameters, especially their dependence on energy res- 
olution, first a and then b in eq. (13). All our results are obtained 
by assuming a reactor of 20 GWth thermal power, a far detector 
of 5 kt fiducial volume with 12% weight fraction of free proton 
and 5 years exposure time. 

First we show the expected energy distributions of the reactor 
antineutrinos in Fig. 1. There are four sets of curves for the 
different baseline lengths, 30, 40, 50 and 60 km, from the top to 
the bottom panel. In each panel, the blue and red curves show 



X^tat ^ J 



A'; 



A' 



NH(IH) 



NH(IH) 



(18) ^Recent study [20] argues that the sensitivity estimation based on Gussian 

distribution can be overestimated, but we wiH not pursue this possibility in this 
study. 
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Figure 1: The energy distributions of reactor antineutrino events after 
20 GW(i,-5kt(12% free proton weight fraction)-5yrs exposure at the baseline 
lengths L = 30, 40, 50 and 60 km, in the top-down order. The blue curves are 
for NH, while the red ones for IH. The red an'ows indicate the energies at which 
the difference due to mass hierarchy vanishes. 



the distributions for NH and IH, respectively. The red arrow 
in each panel shows the antineutrino energy at which the mass 
hierarchy dependent term, the last term in eq. (6), vanishes with 
n = 1 in eq. (8b). Most of the reactor antineutrino events are 
expected to populate the energy range between 1.8 MeV and 
8 MeV. We note here that the difference between the NH and 
IH oscillations is due to the difference of the phase A31 defined 
in eq.(5), as shown in eq.(6). This relative phase difference is 
reversed across the arrowed degeneracy point, as most clearly 
seen in the L = 60 km case. 
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Figure 2: The energy distribution of reactor antineutrinos with baseline length 
L = 30 km and 20 GWth-5kt(12% free proton weight fraction)-5yrs exposure. 
Upper: The case with exact Ey measurement where the dashed blue and dashed 
red curves are for NH and IH, respectively. The solid curve shows the best fit 
of IH assumption to the NH data. The red an'ow points out the energy at which 
the dift'erence due to mass hierarchy vanishes. Lower: 6/ VEw7% energy res- 
olution case. 



Figures 2 and 3 show energy distributions for L = 30 km 
and 50 km, respectively, in which the exact Ey measurement is 
assumed for the upper panel, whereas in the lower panel the en- 
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Figure 3: Same as Fig. 2 but with baseline length L = 50 km 



ergy resolution of a = 6% with b — in eq. (13) is assumed. 
The dashed blue curve corresponds to the NH case, and the 
dashed red curve to the IH case, while the solid curve is ob- 
tained using the parameter values fitted to the NH data with the 
"wrong" IH assumption. At L = 30 km, the solid curve almost 
coincides with the dashed blue one even with the exact energy 
measurement, implying that it is almost impossible to distin- 
guish the mass hierarchy by experiments at L = 30 km. This 
is because the small phase shift between the NH and IH predic- 
tions can be absorbed by a small shift in lAmjjl by a fraction 
of its present uncertainty, 0.1 x lO^^eV^. The situation only 
becomes worse with introducing a finite energy resolution. 

The situation changes when the second peak, the n - 2 point 
in eq. (8a), of the mass hierarchy dependent term appears in the 
energy range. The mass hierarchy difference can no longer be 
absorbed by a shift in |An73j| since the relative phase difference 
between the NH and IH oscillations changes across the degen- 
eracy point. There is no way to make the differences on the both 
sides compensated, resulting in the distinct mismatch between 
the dashed blue curve (for the NH data) and the solid curve 
(the best-fit under the IH assumption) as shown in the upper 
panel of Fig. 3, where the antineutrino energy is exactly mea- 
sured. Once the finite energy resolution is introduced, the phase 
difference in the lower energy side of the degeneracy point is 
significantly smeared out as it oscillates faster w.r.t. Ey at the 
low energy, hence it is easier for one oscillation period to be 
covered by a sizable Gaussian profile of the detector response 
function. The remaining difference in the higher energy side 
can then be absorbed by a small shift in lAmjJ, resulting in an 
excellent fit (solid curve) to the NH data (blue dashed curve) 
in the lower panel of Fig. 3, shown for 6%/ yZs/MeV energy 
resolution. From these result, we can conclude that the physics 
potential for mass hierarchy discrimination strongly depend on 
the energy resolution. 

To discuss more qualitatively the parameter shifts which have 
resulted in the excellent fits, we plot in Fig. 4 the pull factors 
of the five fitting parameters, sin^20i2, sin^20i3, Amjj, lAm^J 
and fsys, as functions of the baseline length L. The pull factor 
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Figure 4: The best-fit values for sin 26i2, sin" 26*13, Am^j, |Amj|| and fsyj, 
v.s. the baseline length L with a = 2%, 3%, 6% and h = 0. The results 
for 20 GWth-5kt(12% free proton weight fraction)-5yrs exposure are shown by 
solid curves for NH, and dashed curves for IH. 



of parameter Y is defined as (y^" - Y"'P"'^ /SY, and its square 
contributes to the function of eq.(15). The best fit values 
with the wrong hierarchy assumption are shown by green, blue 
and red curves for a - 2,3 and 6% with b - 0, respectively. 
As expected, |Am^j| shifts significantly with a negative (NH) 
or positive (IH) pull factor of 0.5 or less, especially in short 
baseline lengths. Although sin^ 20i3 also seems to contribute 
significantly at L ~ 30- 80 km for the a - 2% and 3% cases, we 
checked that its contribution for reducing (Ax^)mm is negligible 
compared to |Amjj|. The other parameters do not contribute 
significantly. At large baseline length, L > 80 km, none of the 
model parameters gives a significant pull factor. 

Figure 5 shows the resulted {Ax^)mw value as a function 
of the baseline length L, for several energy resolutions, a = 
2, 3, 4, 5 and 6% (with b -0)m eq.(13), from the top to the bot- 
tom. Solid curves are for NH, while dashed curves are for IH. 
The results clearly show that the mass hierarchy can be deter- 
mined by those experiments only if the energy resolution of the 
detector is 3% / y/E/MeV or better, and that the optimal base- 
line length (as shown by the cross symbol) is around 50 km 
for that resolution. The small (A;t'^)„„„ for the baseline length 
L < 40 km and L > 80 km is due to a shift in lAm^J and low 
statistics, respectively. For the a - 5 and 6% cases (Ax^)mm 
stays almost zero at all L. 

Next we discuss the effect of the systematic uncertainty part 
of the energy resolution, b, in eq.(13). The Fig. 6 shows the 
{Ax^)mm value as a function of the baseline length L for different 
b values with a-l>%. The curves from the top to the bottom are 
obtained for b = 0%, 0.5%, 0.75% and 1%, respectively. The 
effect of the systematic uncertainty is significant as discussed in 
Ref. [18], reducing the peak value of {Ax^)min from 4.5 {b - 0) 
to 4 (^7 = 0.5%), 3.5 {b = 0.75%) and 3 (fe = 1%) for NH. The 
optimal L shortens from 52 km for {a,b) - (3,0)% to 47 km 
for (fl, b) - (3, 1)%. The Fig. 7 shows another similar figure for 




Figure 5: {fiiX )mm for mass hierarchy discrimination shown as a function of the 
baseline length L, when the energy resolution in eq. (13) is varied with a = 2 
to 6% and h = 0, from the top to the bottom. The results for 20GW,h-5kt(12% 
free proton weight fraction)-5yrs exposure are represented by solid curves for 
NH, and by dashed curves for IH. The cross symbols mark the optimal basehne 
lengths. 
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Figure 6: (^X^)mm for mass hierarchy discrimination v.s. baseline length 
L, with the energy resolution in eq. (13) being a = 3% and b = 
0%, 0,5%, 0.75%, 1%, from the top to the bottom. The other conditions are 
the same as Fig. 5. 
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a - 2%. In this case {hx^)mm is reduced from 13.3 {b - 0) to 



(a, b) = (2, 0) NH 
IH 
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Figure 7: Same as Fig. 6 but with the energy resolution a 
Wo, 0.5%, 0.75%, 1% from the top to the bottom. 



2% and h 



11.8 ib = 0.5%), 10.3 (b = 0.75%) and 8.5 {b = 1%). 

In addition, the neutrino parameters, 

sin^ 2012, sin^ 2013, A/Wjj and lAmjJ, can be measured ac- 
curately with statistical uncertainties shown in Fig. 8. We 
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Figure 8: The statistical uncertainties of the neutrino model parameters 
measured by this experiment as functions of the baseline length L after 
20GWti|-5kt(12% free proton weight fraction)-5yrs exposure. The results for 
both hierarchy (NH by solid and IH by dashed curves) and for the energy reso- 
lution of eq.(13) with (a,b) = (3,0.5), (3, 1) and (6, 1)% are shown. 



find 



(5sin2 26ii2 ~4x 10^^(0.5%), 



21 -3 X 10"^eV^(0.4%), 

31 1 



6Am: 



6\AmU ~7x 10"''eV2(0.3%) 



(21a) 
(21b) 
(21c) 



with the energy resolution of {a,b) = (3,0.5)% at L = 50 
km; the percentage values in the parentheses denote the rel- 
ative accuracy of the measurement. Those uncertainties are 
almost independent of the mass hierarchy and of the energy 
resolution, with the only exception of the lAw^jl uncertainty 
for which the larger resolution results in the larger uncertainty: 



lAm^J ~ 8 X lO^^eV- for the resolution ia,b) = (3, 1)% and 
1.8 X lO^^eV^ for (a, b) = (6, 1)% at L = 50 km. The uncertain- 
ties of sin^ 2012 and Am|j show the rapid reduction after L = 20 
km and stabilize for L > 40 km. This is because the normal- 
ization and shape of the slowly varying oscillation pattern in 
Fig. 1 determine sin^20i2 and Awjj, respectively, which is al- 
most independent of the energy resolution. On the other hand, 
sin^ 2013 and \ Am\^ \ are measured most accurately around L ~ \ 
km, which motivated the first round of the reactor antineutrino 
oscillation experiments such as Daya Bay [1], RENO [2] and 
Double Chooz [25]. The uncertainty of sin^ 20i3 quickly grows 
to the Daya Bay expectation of 5% [29], which is implemented 
as the input in this analysis, at L > 30 km. Somewhat sur- 
prisingly, the uncertainty of lAm^J remains small at the level 
of 1 X lO^^eV^ up to L ~ 60 km when energy resolution is 
3%/ V-E/MeV or better. We find that this is because the rapid 
oscillation pattern due to lAm^J can be resolved even after the 
smearing in the observed energy as can be seen in Fig. 3. With 
better energy resolution, more oscillation patterns are recog- 
nized and higher accuracy of the lAm^J measurement can be 
achieved. 



5. Discussions and Conclusion 

In this paper we have investigated the sensitivity of medium 
baseline reactor electron-antineutrino oscillation experiments 
for determining the neutrino mass hierarchy by performing 
the standard analysis. We carefully study the impacts of 
the systematic uncertainty of the energy resolution {6E/E)^ = 

{a I y/E/MeW^ + b^ and find that the optimal baseline length 
and the sensitivity strongly depend on the energy resolution. 
The optimal baseline length which maximizes the mass hi- 
erarchy resolving power of the experiment is found to de- 
pend slightly on the energy resolution, preferring the length 
slightly shorter than 50 km for the energy resolution of 
(a,b) = (3,0.75)%,(3,1)%,(2,0.75)% and (2, 1)%. At the 
optimal baseline length, the energy resolution better than the 
3%/ V-E/MeV level is needed to determine the neutrino mass 
hierarchy pattern. Three-sigma determination of the mass hier- 
archy is possible for an experiment with 20 GWth-5kt(12% free 
proton weight fraction)-5yrs exposure if an energy resolution 
of (fl, b) - (2, 0.75)% is achieved, while a factor of 3 larger or 
longer experiment is needed to achieve the same goal for the 
energy resolution of {a,b) - (3,0.75)%. It is also found that 
this experiment can measure the neutrino parameters, sin^ 20i2, 
Amjj and |Am^j |, very accurately as shown in (21) for an exper- 
iment of 20 GWth-5kt(12% free proton weight fraction)-5yrs at 
L~50km. 
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